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1.  INTRODUCTION 


During  the  course  of  the  erosion  of  a  gun  tube  bore  surface  by  the 
combustion  products  of  burning  propellants,  significant  alterations  to 
that  surface  are  caused  by  the  interaction  between  the  metal  and  the 
hot  gases  at  high  pressures.  After  only  a  small  number  of  exposures, 
physical  and  metallurgical  changes  to  the  surface  are  manifested  by 
the  appearance  of  cracking  along  the  grain  boundaries  and  a  heat 

affected  region  known  as  the  "white"  layer^.  The  surface  most  likely 
undergoes  the  austentite/martensite  transformation  with  each  exposure 
causing  volume  changes  which  result  in  heat  checking  and  deep  cracking 

after  many  exposures'*".  In  many  cases  flow  patterns  are  observed,  in¬ 
dicating  that  the  surface  melts  and  then  is  sheared  by  the  flowing 
2  3 

gases  *  .  In  addition  to  the  metallurgical  and  physical  changes  taking 
place,  one  can  also  expect  to  find  chemical  effects  which  are  present 
only  in  the  surface  layers.  These  effects  are  due  to  chemical  reactions 
between  the  combustion  products  and  the  metal  at  the  surface  and/or 
physical  effects  such  as  enhanced  diffusion  along  grain  boundaries.  Al¬ 
though  gross  changes  have  been  well  documented  through  the  use  of  instru¬ 
ments  such  as  Scanning  Electron  Microscopes  (SEM) ,  borescopes,  hardness 
testers,  and  a  variety  of  other  metallurgical  methods,  the  instrumentation 
and  techniques  capable  of  measuring  the  chemical  changes  on  the  surface 
have  only  been  applied  recently.  Techniques  such  as  X-ray  Photoelectron 
Spectroscopy  (XPS) ,  Characteristic  X-ray  Analyses  and  Ion  Beam  Analysis 
(IBA)  are  now  being  applied  to  the  study  of  gun  bore  surfaces. 

To  date,  most  experiments  on  the  effects  of  erosion  on  gun  steel 
have  been  done  in  tests  where  large  numbers  of  rounds  have  been  fired 
through  a  barrel.  Afterwards,  the  bore  surface  was  destructively  examined. 

In  such  experiments,  gross  effects  such  as  heat  checking,  fatigue  cracking 
and  in  cases  where  additives  are  used,  barrel  fouling  may  obliterate  many 
of  the  details  of  the  erosion  process.  A  method  by  which  the  outer  sur¬ 
face  layers  of  a  steel  sample  exposed  to  the  erosive  effects  of  burning 
propellants  can  be  examined  between  consecutive  shots  would  be  very  use¬ 
ful.  Since  in  a  typical  large  caliber  gun  the  erosive  wear  amounts  to 
about  1  ym  per  shot,  it  is  highly  probable  that 

1.  Sidney  Breitbart,  "On  the  Erosion  of  Metallic  Nozzles  by  Powder 
Gases,"  Ballistic  Research  Laboratory  Memorandum  Report  No.  533, 

April  1951.  AD  #802148 

2.  Franklin  A.  Vassalo  and  W.  Richard  Brown,  "Shock  Tube  Gun  Melting 
Erosion  Study,"  Contract  Report,  ARBRL-CR-00406,  August  1979.  AD  #A076219 

3.  A.C.  Alkidas,  S.O.  Morris,  C.W.  Christoe,  L.H.  Caveny,  and  M.  Sum- 
merfield,  "Erosive  Effects  of  Various  Pure  and  Combustion-Generated 
Gases  on  Metals,  Part  II,"  Contract  Report,  AMMRC  CTR  77-25, 

October  1977. 
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such  individual  shots  affect  a  surface  layer  of  approximately  the  same 
thickness.  Thus,  this  surface  analysis  method  must  be  capable  of  exam¬ 
ining  that  outer  1  ym  layer  of  the  steel. 

At  the  BRL,  a  surface  analysis  technique  employing  ion  beams  has  been 
developed  which  is  capable  of  quantitative  measurement  of  the  concentra¬ 
tions  and  depth  profiles  of  carbon,  nitrogen  and  oxygen  in  the  outer  mi¬ 
cron  layers  of  eroded  samples.  These  elements  are  of  greatest  interest 
because  their  compounds  constitute  a  large  percentage  of  a  given  propel¬ 
lant.  As  the  propellant  combustion  products  interact  with  the  steel  sur¬ 
face,  they  either  melt  the  surface  and  blow  it  away  or  change  the  surface 
layer  compositions.  When  additives  are  used,  interactions  may  also  occur 
between  the  hot  propellant  gases  at  high  temperatures  and  the  additive 
compounds.  This  raises  the  possibility  that  the  heated  surface  acts 
either  as  a  catalyst  or  condenser  of  the  gaseous  products.  A  measure  of 
the  carbon,  nitrogen  and  oxygen  concentrations  left  behind  and  the  depth 
profile  of  the  affected  layer  in  the  steel  surface  can  be  expected  to 
yield  detailed  information  about  the  nature  of  the  processes  involved  be¬ 
tween  the  hot  gases  and  the  surface  during  the  propellant  burn. 

A  number  of  experiments1,4,5  have  indicated  that  erosion  reduction 
may,  in  part,  be  due  to  the  formation  of  protective  layers  on  the  steel 
surface  which  are  often  quite  thin.  These  layers,  of  different  physical 
and  chemical  composition  than  the  substrate  steel,  protect  that  steel  by 

reducing  the  flow  of  heat  to  the  steel6.  The  characteristics  of  such 
protective  layers,  the  conditions  under  which  they  form  and  other  details 
needed  for  their  optimization  are  not  understood.  In  addition,  although 
the  erosion  reducing  properties  of  propellant  additives  such  as  TiO^  wax 

have  long  been  appreciated,  the  detailed  mechanisms  by  which  these  addi¬ 
tives  work  have  not  been  clearly  described.  For  a  complete  review  of 
efforts  to  reduce  gun  tube  erosion  by  the  use  of  high  temperature  metal 
coatings  as  well  as  various  wear  reducing  additives,  see  the  Proceedings 

7 

of  the  Tri-Service  Gun  Tube  Wear  and  Erosion  Symposium. 

4"!  J.  Richard  Ward  and  Timothy  L.  Brosseau,  "Role  of  the  Insulating 
Layer  From  TiC^-Wax  Liner  in  Reducing  Gun  Tube  Wear,"  Ballistic 

Research  Laboratory  Technical  Report  ARBRL-TR-02238,  April  1980. 

5.  R.  Birkmire  and  A.  Niiler,  "Applications  of  the  Radioisotope  Wear 
Measurement  Technique,"  Ballistic  Research  Laboratory  Technical 
Report  ARBRL-TR-02075 ,  June  1978.  AD  #A058307 

6.  L.H.  Russell,  "Simplified  Analysis  of  the  Bore  Surface  Heat  Transfer 
Reduction  in  Gun  Barrels  as  Achieved  by  Using  Wear-Reducing  Additives," 
Naval  Surface  Weapons  Center  Report  TR-3378,  October  1975. 

7.  "Proceedings  of  the  Tri-Service  Gun  Tube  Wear  and  Erosion  Symposium," 
edited  by  Jean-Paul  Picard  and  Iqbal  Ahmad,  March  1977. 
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The  work  that  will  be  described  in  this  report  was  undertaken  in 
order  to  develop  data  on  the  effects  that  the  burning  of  various  pro¬ 
pellants  in  combination  with  wear  reducing  additives  have  on  the  surface 
layer  compositions  of  a  typical  gun  steel,  AISI  4340.  Central  to  this 
study  was  the  realization  that  ion  beam  analysis  methods  could  be  used 
to  determine  the  properties  of  bore  surfaces  that  bear  most  strongly  on 
^  erosion.  The  experiments  were  designed  so  that  shot  by  shot  trends  and 

correlations  could  be  observed  in  a  small  number  of  rounds.  Correlating 
these  trends  to  existing  wear  data  on  widely  used  propellant/additive 
systems,  leads  to  predictions  concerning  new  and  novel  systems  on  the 
basis  of  relatively  small  expenditures  of  time,  effort  and  money.  These 
experiments  had  to  be  carried  out  in  the  laboratory  where  the  maximum 
number  of  parameters  could  be  controlled. 


2.  EXPERIMENTAL  PROCEDURE 
2 . 1  Test  Firing  Program 


In  order  to,  control  the  largest  possible  number  of  experimental 
parameters  as  well  as  to  allow  non-destructive  analysis  of  the  steel 
surfaces  exposed  to  burning  propellants,  contoured  nozzles  of  4340  steel 
were  used.  The  4340  steel,  hardened  to  42  Rockwell,  closely  resembles 
the  composition  and  hardness  of  typical  gun  tube  steel.  Twelve  of  these 
nozzles  were  exposed  in  a  37  mm  gun  whose  barrel  was  cut  off  and  which  was 
modified  to  act  as  a  blow-out  chamber.  A  sketch  of  this  37  mm  chamber 
is  shown  in  Figure  1.  The  1/16"  thick  mild  steel  discs  were  used  to 
contain  the  burning  propellant  gases  until  the  pressure  became  high 
enough  to  rupture  the  disc  and  send  the  hot  gases  streaming  over  the  noz¬ 
zle  surface  causing  the  erosion.  Each  of  twelve  nozzles  was  fired  with 
a  different  combination  of  the  propellants  M2,  M30  and  Ml  and  the  addi¬ 
tives  TiO^  wax  (XM1) ,  Talc  Wax  (TW) ,  Polyurethane  Foam  (PUF)  and  no  addi¬ 
tive  (NA) .  The  main  criterion  for  choosing  these  propellants  was  the  large 
range  of  their  flame  temperatures  from  2417°K  for  Ml,  3040°K  for  M30  to 
3320°K  for  M2. 

With  two  exceptions,  four  shots  were  fired  through  each  nozzle.  The 
propellant  charge  was  chosen  to  give  approximately  the  same  pressure-time 
trace  with  each  of  the  three  propellants.  The  pressures  were  measured  in 
the  chamber  with  a  Top-Hat  strain  gauge  and  the  pressure-time  trace  was  re¬ 
corded  for  each  shot  on  film.  A  sample  pressure-time  trace  is  shown  in 

g 

Figure  2.  The  M2  charge  mass  was  taken  directly  from  previous  BRL  work 
and  the  others  were  matched  to  it  by  the  pressure-time  trace.  The  amount 
of  the  additives  which  were  used  was  also  taken  directly  from  ref.  8. 

All  the  additives  came  in  sheet  form,  were  rolled  around  the  forward  part 

8.  Michael  A.  Schroeder  and  Masahiro  Inatome,  "The  Relationship  Be¬ 
tween  Chemical  Composition  and  Wear-Reducing  Effectiveness  of  Some 
Laminar  Additives  for  Gun  Propellants:  Polyvinyltetrazole," 

Ballistic  Research  Laboratory  Memorandum  Report  2512,  August  1975. 

AD  #B007029L 
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1.  RETAINING  RING 

2.  NOZZLE 

3.  BLOWOUT  DISC 

4.  SPACER 

5.  O-RING  SEAL 


Figure  1.  37mm  Vented  Chamber. 
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Figure  2.  Pressure-time  trace  for  nozzle  #10. 
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of  the  propellant  bed  and  the  ends  were  flapped.  This  configuration  of 
the  additive  sheet  has  been  shown  to  be  the  most  effective  at  erosion  re- 

duction  '  .  The  propellant/additive  parameters  under  which  all  of  the 
twelve  nozzles  were  fired  are  shown  in  Table  I.  The  additive  masses 
shown  are  averages  for  all  shots  taken  with  the  given  configuration. 

After  each  shot,  the  nozzles  were  cleaned  with  a  water-detergent 
wash,  rinsed  in  alcohol  and  the  mass  losses  measured.  The  record  of  the 
pressure  and  mass  loss  measurements  is  shown  in  Table  II.  The  averages 
and  their  mean  standard  deviations  are  shown  for  the  peak  pressures,  the 
durations  of  the  pressure  pulses  and  the  mass  losses  for  each  of  the  pro¬ 
pellant/additive  combinations.  The  variation  in  the  peak  pressure  was 
found  to  be  up  to  10%  and  is  most  likely  due  to  the  variation  in  the  man¬ 
ner  in  which  the  steel  blowout  discs  ruptured.  The  small,  but  clear, 
increase  in  the  pressures  when  additives  were  used  was  probably  due  to  the 
increase  in  the  amount  of  gas  from  the  burning  additives  as  well  as  the 
lowering  of  the  effective  volume  of  the  chamber.  The  duration  of  the  pro¬ 
pellant  burn  seems  to  be  mainly  a  function  of  the  propellant  type.  It  is 
significant  to  the  extent  that  a  longer  time  at  a  temperature  above  sur¬ 
face  melting  point  produces  more  erosion.  However,  since  no  attempt  was 
made  to  measure  the  nozzle  surface  temperatures  in  this  experiment,  the 
pulse  duration  does  not  contribute  significantly  in  our  conclusions. 

The  shot-to-shot  variation  in  the  mass  losses  indicated  by  the  stand¬ 
ard  deviations  in  Table  II,  was  somewhat  larger  than  the  measurement  un¬ 
certainties  of  ±  0.1  mg  would  allow.  However,  other  experiments  with  this 

3 

chamber  as  well  as  with  other  test  fixtures  have  shown  very  similar  shot- 
to-shot  variations.  Although  the  reason  for  this  variation  is  not  known, 
we  speculate  that  localized  random  variations  in  the  flame  patterns  can 
be  responsible.  Another  possibility  is  that  the  blow-out  discs  may  have 
sheared  unevenly,  providing  uneven  regions  of  hot  gas  flow. 

2 . 2  Ion  Beam  Analysis 


After  the  cleaning  and  mass  loss  measurements,  the  nozzles  were 
subjected  to  ion  beam  analysis.  Most  of  the  details  of  the  technique 

are  fully  described  in  an  earlier  report^  and  so  will  not  be  discussed. 
However,  since  a  different  accelerator  was  used  in  these  studies,  operat¬ 
ing  at  a  different  beam  energy,  there  are  some  differences  in  details 
which  will  be  pointed  out  when  appropriate.  It  is  recommended  that  in 
order  to  fully  understand  the  IBA  method,  Ref.  10  and  references  therein 
should  be  studied. 


9^  Timothy  L.  Brosseau  and  J.  Richard  Ward,  "Reduction  of  Heat  Transfer 
in  105  mm  Tank  Gun  by  Wear-Reducing  Additives,"  Ballistic  Research 
Laboratory  Memorandum  Report,  BRL-MR-2698,  November  1976.  AD  #B015308L 
10.  A.  Niiler,  J.E.  Youngblood,  S.E.  Caldwell  and  T.J.  Rock,  "An 

Accelerator  Technique  for  the  Study  of  Ballistic  Surfaces,"  Bal¬ 
listic  Research  Laboratory  Report  No.  1815,  August  1975.  AD  #A016899 
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Table  1.  Firing  Data 


Nozzle 

Numb. /Add./ 
Rounds  Fired 

Propel  1  ant 
Type/Mass/Temp. 

Peak 

Pressure 

(MPa) 

Mass  Loss 
(mg) 

2/NA  /4 

165. 5±  3.4 

36. 8i5. 5 

3/XM1/4 

M2/85g 

180. 0±1 1.0 

8. 5±1 . 1 

4/TW  /4 

3320° K 

1  78. 6±1 3. 1 

6. 9i3. 2 

5/PUF/4 

174. 5±  3.4 

9. 9i5. 0 

6/NA  /4 

1  78.  6±  5.5 

1 3. Oil .  0 

7/XM1/4 

M30/1 1 Og 

186. 2±  4.8 

5.9i3.6 

8/TW  /5 

3040° K 

182.1 ±17. 2 

4.6il .9 

9/PUF/4 

187. 6±  4.8 

6. 5i2. 5 

10/NA  /4 

169. 0±  3.4 

4.9il .9 

1 1/XM1/3 

Ml/1 05g 

1 79. 3±1 5. 2 

4. Oil . 1 

12/TW  /4 

■  241 7° K 

180.  Oil  4. 5 

3. 9il . 2 

1 3/PUF/4 

179. 3i  4.8 

3.4il . 9 

Table  2.  Results 

of  Pressure  and 

Mass  Loss  Measurements 

Nozzle 

Propel  1  ant/ 

Peak 

Duration 

Mass  Loss 

Additive 

Pressure 

(K  psi) 

FWHM  (msec) 

(mg) 

2 

M2/NA 

24. OiO. 5 

2. 20i0. 1 0 

36.8i5.5 

3 

M2/XM1 

26. 1 il . 6 

2. 08i0. 16 

8. 5+1.1 

4 

M2/TW 

25. 9il . 9 

2. 06i0. 1 6 

6.9i3.2 

5 

M2/PUF 

25. 3i0. 5 

2. 14±0. 10 

9. 9i5. 0 

6 

M30/NA 

25. 9i0. 8 

2. 54i0. 04 

13. Oil. 0 

7 

M30/XM1 

27. OiO. 7 

2. 48i0. 04 

5. 9i3. 6 

8 

M30/TW 

26. 4i2. 5 

2. 68i0. 24 

4.6il .9 

9 

M30/PUF 

27. 2i0. 7 

2. 50i0. 08 

6. 5i2. 5 

10 

Ml/NA 

24. 5i0. 5 

2. 36i0. 06 

4.9il .9 

11 

M1/XM1 

26. 0i2. 2 

2. 26i0. 14 

4. Oil  .1 

12 

Ml/TW 

26. Ii2. 1 

2. 28i0. 1 2 

3.9il .2 

13 

Ml/PUF 

26. OiO. 7 

2. 28i0. 1 0 

3.4il  .9 

The  IBA  method  is  outlined  as  follows:  A  high  velocity,  focused 
beam  of  ions,  deuterons  in  this  case,  is  generated  by  an  accelerator, 
and  directed  onto  the  surface  of  a  sample  to  be  analyzed.  As  the  deuter¬ 
ons  impinge  on  the  sample,  they  can  undergo  nuclear  reactions  as  well  as 
elastic  scattering  with  the  nuclei  of  the  sample1 s  constituents.  The 
type  of  particles  that  emerge  from  these  encounters,  as  well  as  the  ener¬ 
gies  of  these  particles,  is  very  specific  to  the  type  of  target  nucleus  i 

with  which  the  encounter  takes  place.  In  addition,  the  number  of  a  given 
type  of  emitted  particles  as  well  as  its  observed  energy  distribution  can 
be  used  to  deduce  concentrations  and  depth  profiles  of  the  nuclei  in  the 
surface  with  which  the  interactions  took  place.  Figure  3  is  used  to  illus¬ 
trate  this  method.  The  insert  in  this  figure  shows  a  surface  containing 
iron,  carbon,  nitrogen  and  oxygen.  As  the  deuteron  beam  enters  the  surface, 
one  of  the  many  different  types  of  encounters  the  deuteron  particles  under¬ 
go  is  elastic,  billiard  ball  type,  scattering  from  any  of  the  constituents. 

If  the  scattered  particles  are  detected  at  an  angle  greater  than  90°  with 
respect  to  the  beam  direction,  this  elastic  interaction  is  known  as  Ruther¬ 
ford  Backscattering  (RBS) .  The  spectrum  shown  in  Figure  3  shows  this  RBS 
distribution  to  the  left.  Other  types  of  encounters  between  the  beam  deu¬ 
terons  and  target*  are  exothermic  nuclear  reactions  (NR)  with  the  carbon, 
nitrogen  and  oxygen  in  the  surface.  In  these  reactions,  higher  energy 
protons  and  a-particles  are  emitted.  The  proton  peaks  due  to  C,N  and  0  are 
also  shown  in  figure  3  at  the  right.  In  Ref.  10,  the  NR  analysis  was  done 
on  a-particles  emitted  from  oxygen. 

Both  the  RBS  distribution  and  the  NR  peaks  can  be  used  to  obtain 
specific  information  about  the  surface  layers.  The  NR  peaks  are  very 
specific  in  identifying  the  constituents  of  the  surface  while  the  shape 
of  the  RBS  distribution  is  very  sensitive  to  the  total  amount  of  ’’contamina¬ 
tion”  on  the  surface.  The  most  detailed  information  about  concentra¬ 
tions  and  depth  profiles  of  C,N  and  0  is  obtained  from  the  NR  peaks. 

The  area  under  a  given  peak  is  related  to  the  concentration  of  that 

element  in  the  surface  while  the  shape  of  a  peak,  especially  its 

width,  is  related  to  the  depth  profile.  In  order  to  quantitatively 
determine  these  concentrations  and  depth  profiles,  a  computer  program 
which  models  the  experimental  spectrum  by  numerical  integration,  is 
used.  This  program  will  be  described  in  a  later  section. 

A  deuteron  beam  of  nominal  energy  of  0.665  MeV,  obtained  from  the 

BRL  Cockcroft-Walton  accelerator  was  used  for  these  analyses.  The  beam 

energy  was  checked  periodically  and  referenced  to  the  0.669  MeV  resonance  * 

in  the  19F(p,ay)160  reaction.  The  uncertainty  in  the  beam  energy  was  found 
to  be  ±  .003  MeV.  Samples  were  mounted  in  a  scattering  chamber  normally 

-  6  ^ 
held  at  a  pressure  of  2  X  10”  torr.  A  silicon  surface  barrier  detector 

was  used  for  the  simultaneous  detection  of  both  the  RBS  and  NR  products. 

This  type  of  detector  produces  an  output  current  pulse  which  is  propor¬ 
tional  to  the  energy  of  the  incident  charged  particle.  These  current  pul¬ 
ses  are  passed  through  a  charge  sensitive  preamplifier,  a  high  resolution 
linear  amplifier  and  finally  stored  in  a  multi-channel  pulse  height  analyzer. 
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DEUTERQN  BEAM 
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Figure  3.  Ion  Beam  Analysis  Schematic  and  Spectrum. 


The  RBS  distribution  and  NR  peaks  shown  in  figure  3  result  from  this  type 
of  signal  analysis.  In  order  to  obtain  the  maximum  detector  resolution 
of  0.012  MeV  for  5.5  MeV  alphas,  a  pulse  pile-up  rejection  feature  is 
incorporated  into  the  electronics.  In  a  pulse  counting  experiment  such  as 
this,  it  is  possible  for  two  independently  generated  particles  to  produce 
overlapping  pulses  in  the  detector/preamplif ier/amplif ier/analyzer  chain. 
When  such  a  double  event  occurs,  the  total  system  resolution  is  degraded 
and  under  some  conditions,  the  spectral  shapes  can  also  be  distorted.  The 
pulse  pile-up  rejection  circuitry  used  eliminated  about  95-6  of  these  double 
events  as  long  as  the  two  pulses  were  more  than  75  nsec  apart  and  produced 
by  particles  with  energies  greater  than  50  keV.  The  samples  were  enclosed 
by  an  electron  suppression  shield  permitting  accurate  beam  charge  measure¬ 
ments.  The  beam  current  integrator  used  to  measure  beam  flux  was  gated 
to  correct  for  the  electronic  system  dead  time. 


3.  ANALYSIS  PROCEDURE 


3 . 1  Program  PROFILE 

The  program  PROFILE  was  used  to  analyze  the  RBS  distributions  and 
the  NR  peaks  from  the  nozzles  in  order  to  obtain  the  concentrations 
and  depth  profiles  of  the  surface  contaminants.  A  complete  description 
of  the  program  as  it  was  originally  used  is  given  in  Ref.  10  but  since 
significant  changes  have  been  incorporated  since  then,  a  detailed 
description  will  also  be  given  now.  Additional  information  including 
a  listing  and  a  cataloging  of  the  description  for  its  use  on  the  BRL 
CDC  7600  computer  is  given  in  reference  11. 

The  procedure  used  in  the  calculation  is  essentially  an  iterative 
fitting  of  the  experimentally  measured  spectrum.  The  sample  surface  is 
divided  into  a  number  of  layers  (up  to  100),  each  of  which  is  assigned  a 
thickness  and  a  specific  composition  which  may  include  up  to  five  types 
of  atoms.  All  of  the  experimental  parameters  such  as  beam  energy  and 
intensity,  detector  geometry  and  angle  are  known,  while  the  specific  com¬ 
positions  of  the  layers  are  the  only  unknowns.  These  latter  quantities 
are  varied  in  proceeding  to  an  acceptable  fit  to  the  data.  Although  in 
the  final  analysis  the  layer  compositions  are  the  basic  unknowns  in  the 
experiment,  there  are  a  number  of  prior  clues  which  eliminate  the  need  for 
a  blind  search  through  the  whole  periodic  table.  First,  there  is  some  prior 
knowledge  of  the  surface.  The  4340  steel  is  over  95%  iron  and  contains  up 
to  .4%  carbon.  The  fact  that  the  surface  was  exposed  to  propellants 
would  also  indicate  that  the  propellant  constituents  such  as  carbon, . 
nitrogen  and  oxygen  as  well  as  some  additives  such  as  sulfur  and  barium 
from  BaS04  might  be  expected.  Second,  the  presence  of  certain  peaks  in 

IT!  A.  Niiler,  R.  Birkmire  and  J.  Gerrits,  "PROFILE:  A  General  Code 

for  Fitting  Ion  Beam  Analysis  Spectra,"  Ballistic  Research  Laboratory 
Technical  Report  ARBRL-TR-02233,  April  1980.  AD  #A084984 
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the  NR  region  are  definitive  signatures  of  the  presence  of  lighter  elements 
such  as  carbon,  nitrogen  and  oxygen.  Finally,  if  there  is  enough  of  a  given 
light  element  on  the  surface,  it  will  produce  a  peak  on  top  of  the  RBS  dis¬ 
tribution  from  iron. 

The  calculation  proceeds  by  using  layer  constituents  which  result 
from  "educated"  guessing  and  determining  the  specific  energy  loss  tables 
as  a  function  of  the  layer,  type  of  beam  ion  and  the  ion  energy.  The 
reference  data  for  energy  losses  are  obtained  from  the  tabulations  of 

12 

Anderson  and  Ziegler  .  The  accuracy  of  the  reference  energy  losses  is 
very  important  because  as  the  ions  (deuterons,  protons,  ot-particles)  pass 
through  the  layers,  they  lose  energy.  Their  final  detected  energy  distri¬ 
butions  are  very  sensitively  dependent  on  the  values  of  the  energy  losses 
as  obtained  from  the  reference  tables. 

At  this  point,  each  of  the  defined  layers  of  the  surface  is  sub¬ 
divided  into  slabs  in  order  to  improve  the  energy  resolution  of  the 
calculation  and,  in  turn,  each  slab  is  illuminated  by  the  same  unit  flux 
of  deuterons.  The  deuterons  are  then  required  to  undergo  a  particular 
type  of  interaction,  either  RBS  or  NR,  with  one  of  the  layer  consti¬ 
tuents.  The  deuteron  energies  vary  from  one  layer  and  slab  to  the  next 
depending  on  how  much  energy  loss  has  taken  place  since  the  beam  entered 
the  surface.  In  this  calculation,  the  RBS  and  NR  interactions  are  run 
separately,  but  contributions  from  all  target  constituents  are  added 
for  each  type  of  interaction.  The  probability  for  a  given  type  of 
interaction  is  determined  from  the  cross  sections  for  those  reactions. 

In  all  cases,  the  RBS  probabilities  are  calculated  from  the  Rutherford 
elastic  cross  sections  whose  values  are  given  by 


where  E  is  the  deuteron  energy, 

0  is  the  scattering  angle, 

Zj  is  the  deuteron  charge, 

Z^  is  the  target  charge, 
is  the  deuteron  mass, 
and  is  the  target  mass. 

12.  H.H.  Anderson  and  J.F.  Ziegler,  Hydrogen  Stopping  Powers  and  Ranges 
in  All  Elements,  Pergamon  Press,  New  York,  1977. 

13.  N.  Longequeque,  H.  Beaumevieille,  E.  Ligeon,  J.P.  Longequeque,  and 

M.  Sandon,  "Etude  des  Reactions  ^0(d,aQ) , ’^O(d,po)  et  ^0(d,p^) 

de  300  keV  a  1  MeV  (Resultats  Experimentaux) ,  Le  Journal  de  Physique, 
26,  1965,  367. 
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16  17 

The  NR  cross  sections  for  oxygen,  the  0(d,p^)  0  values  were 

13  12  13 

taken  from  Longequeue  et  al.  .  For  the  C(d,Pp)  C  reaction,  we  used 

14  14  15 

values  from  Huez  et  al.  .  For  the  N(d,p<.)  N,  we  used  values  mea¬ 
sured  on  the  BRL  Cockcroft  Walton  accelerator  and  the  University  of 
Delaware  2.5  MeV  Van  de  Graaff.  A  description  of  this  measurement  is 
given  in  Section  4.  The  actual  values  of  these  three  cross  sections 
which  were  used  for  the  oxygen,  nitrogen  and  carbon  analyses  are  shown 
in  figure  4.  The  subscript  on  the  outgoing  particles  (p^,  p^,  p^)  refer 

to  the  excited  unclear  state  of  the  residual  nucleus  in  the  reaction. 

Once  the  deuterons  have  interacted  with  the  target  nuclei  in  a 
given  slab,  the  outgoing  particles  (either  deuterons  in  the  case  of 
RBS  or  protons  in  the  case  of  NR)  are  followed  back  out  of  the  target, 
again  undergoing  the  proper  energy  losses.  Energy  broadenings  due  to 
beam  energy  resolution,  energy  loss  straggling  and  detector  resolution 
are  incorporated  in  the  calculation  at  the  appropriate  points.  When  a 
distribution  from  a  given  slab  has  been  followed  to  its  exit  from  the 
surface,  it  is  weighted  by  the  cross  section  for  the  appropriate  tar¬ 
get  constituent  at  the  proper  energy  in  the  slab  as  well  as  the  concen¬ 
tration  of  that  constituent  in  the  slab.  The  final  energy  distribution 
is  made  up  of  a  series  of  gaussians  whose  centroids  are  determined  by 
energy  losses,  whose  widths  are  determined  by  energy  broadenings  and 
whose  weights  are  determined  by  the  cross  sections  and  layer  concen¬ 
trations.  These  calculated  distributions  are  fit  by  varying  a  scale 
factor  to  the  experimental  spectra  by  least  squares  methods.  The  concen¬ 
tration  profiles  are  then  adjusted  and  the  calculation  rerun  until  a  sat¬ 
isfactory  fit  is  achieved.  The  primary  criteria  for  acceptable  fits  are 
a  general  agreement  between  the  calculated  and  experimental  shapes  along 
with  the  lowest  possible  chi  squared  for  the  full  calculation.  The  chi- 
squared  is  given  by 


where  N  is  the  number  of  points  fit,  I  is  the  number  of  degrees  of 
freedom,  is  the  calculated  distribution,  is  the  experimental  dis¬ 
tribution  and  A  is  the  least  squares  proportionality  constant  between 
the  two  distributions.  Typically,  I  is  given  by 

I  =  N£(Nc  -  1)  (3) 

where  N0  is  the  number  of  layers  and  N  is  the  number  of  constituents 
X/  c 

14 .  M.  Huez,  L.  Quaglia  and  G.  Weber,  "Fonction  D'Excitation  de  la 

Reaction  ^C(d,pQ)^C  Entre  400  et  1350  keV  -  Distribution  Angulaires," 
Nucl.  Inst,  and  Meth.,  105,  1972,  197. 
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Figure  4.  (d,p)  Cross  Sections  Used  In  Analysis. 
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per  layer.  Since  usually  it  is  possible  to  obtain  acceptable  fits 
within  5  to  10  iterations  of  the  profiles  by  hand,  no  automatic  search 
routines  have  been  incorporated.  Computer  costs  would  be  driven  up 
considerably  by  such  auto-fit  routines. 

The  atomic  concentrations  of  the  various  constituents  of  the  sur¬ 
face  are  obtained  in  a  straightforward  manner  from  the  profiles  using 
the  formula 


N„ 


N. (at/cm  )  =  Y  T.  A 
1 K  J  h  io 

3 


n.  . 
1 3 


N 

l  V* 


(4) 


where  N. 

l 


.  th 


is  the  atomic  concentration  for  i  constituent 
of  the  surface, 

T.!s  are  the  layer  thicknesses. 


is  Avogadros  number. 


th 


n^j!s  are  the^jitomic  percentages  of  the  k  constituent 


in  the  j  layer, 


th 


is  the  atomic  mass  of  the  k  constituent  and 


N0  and  N  are  as  previously  defined. 

A/  C 

The  Tj's  are  program  inputs  whose  values  are  chosen  to  yield  the  desired 

energy  resolution  for  the  calculation.  Thus,  the  only  input  set  which 

determines  the  final  concentration  values  is  the  n,  .'s,  the  set  of 

kj 

atomic  percentages  of  the  layers. 

3.2  Discussion  of  Calculation 


To  understand  the  advantages  and  limitations  of  this  calculational 
procedure  in  obtaining  absolute  concentrations  and  depth  profiles,  ex¬ 
periments  done  on  a  sample  other  than  a  nozzle  will  now  be  discussed. 

For  this  test  case  the  sample  is  a  polished  flat  slab  of  steel  which  was 
heated  to  produce  an  oxygenated  surface  layer.  This  sample  was  then 
analyzed  by  a  .667  MeV  deuteron  beam  resulting  in  RBS  and  NR  spectra  simi¬ 
lar  to  those  shown  in  figure  3.  Using  the  PROFILE  analysis  described  in 
the  previous  section,  the  concentrations  for  C,N,0,  and  Fe  were  determined 
and  are  shown  in  figure  5.  The  oxygen  and  carbon  are  surface  peaked 

2 

with  both  of  them  extending  to  the  maximum  probing  depth  of  .55  mg/cm  . 

There  is  very  little  nitrogen  on  this  surface  extending  only  to  the  .15 

2  2 
mg/cm  depth  and  the  iron  reaches  its  bulk  value  at  the  same  .15  mg/cm 

depth.  The  solid  curve  shown  in  figure  6  is  the  fit  to  the  RBS  and  NR 

spectra  obtained  with  the  set  of  concentrations  shown  in  figure  5.  The 

chi-squared  per  point  for  this  fit  was  1.33. 


22 


100-1 


o 


60- 


i 

g 


-  CARBON 

-  OXYGEN 

- NITROGEN 

- IRON 


ANALYSIS 

DEPTH 

LIMIT 


“  I 

0.4  0.5 


T" 

0.6 


DEPTH  (mg/cm2) 


l 

0.7 


Figure  5.  Concentration  Profiles  of  0,N,C  and  Fe  in  Polished,  Heated, 
Steel  Slab. 
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Figure  6.  Fits  to  the  RBS  and  NR  distributions  from  the  polished,  heated  steel  slab  obtained  with 
various  concentration  profile  sets.  See  text  for  explanation  of  the  different  curves. 


To  test  the  validity  of  IBA  and  the  use  of  PROFILE,  several  checks 
have  been  performed  on  this  sample.  The  sensitivity  of  the  fit  to  varia¬ 
tions  in  the  profiles  was  examined.  To  determine  the  accuracy  of  this 
procedure  for  determining  the  concentration  profile  for  a  given  constitu¬ 
ent,  the  carbon  concentrations  were  lowered  by  10%  over  the  full  analyzed 
depth.  The  dashed  curve  in  figure  6  shows  the  result  of  this  calculation 
for  which  the  chi  squared  was  4.21.  The  largest  contribution  to  the  in¬ 
creased  chi  square  came  from  the  fit  to  the  RBS  spectrum  in  the  regions 
around  .33  MeV  and  .56  MeV.  There  was  almost  no  degradation  in  the  fit 
to  the  carbon  (d,p)  peak  in  the  NR  region  while  the  oxygen  (d,p)  peak  was 
significantly  affected.  This  illustrates  that  the  profiles  for  all  con¬ 
stituents  must  be  correct  in  order  for  the  overall  fit  to  be  good.  It 
further  demonstrates  that  the  calculational  accuracy  in  determining  the 
concentrations  is  better  than  10%.  However,  the  calculated  yields,  and 
therefore  the  concentrations,  are  directly  proportional  to  the  specific 
Rutherford  or  Nuclear  cross  sections  and,  thus,  are  sensitive  to  the  un¬ 
certainties  of  these  cross  sections.  The  Rutherford  cross  section  values 
as  computed  from  equation  1,  are  within  2-10%  of  experimental  values  in 
this  energy  range  while  the  measured  NR  values  have  at  least  a  10%  abso¬ 
lute  uncertainty.-  Consequently,  the  cross  sections  contribute  a  10%  abso 
lute  uncertainty  to  the  concentrations.  The  total  uncertainty  on  the  con 
centration  values  is  therefore  somewhere  between  10  and  15%. 

A  second  variation  in  the  profiles  was  examined  to  determine  the 
depth  resolution  capabilities  of  this  technique.  The  dot-dashed  curve 
in  figure  6  is  the  result  of  increasing  the  thickness  of  the  oxygen 

2 

layer  from  .01  to  .02  mg/cm  .  The  chi-squared  for  this  fit  was  5.94 
with  the  poorer  fit  coming  for  the  RBS  distribution  near  .4  and  .56  MeV 
and  also  in  the  region  of  the  oxygen  peak  in  the  NR  spectrum. 

It  has  been  shown^  that  surface  roughness  of  the  order  of  .02  to 
1.0  microns  will  affect  the  RBS  distributions  if  very  small  incident  or 
exit  angles  are  used.  No  effect  was  found  for  incident  or  exit  angles 
above  5°  or  for  scattering  angles  much  larger  than  100°.  As  in  all  of 
our  experimental  spectra  the  scattering  angle  is  always  160°  and  the 
incident  and  exit  angles  are  always  greater  than  20°,  no  effects  due  to 
surface  roughness  should  be  observable  in  our  spectra.  To  illustrate 
this  fact,  we  produced  several  surfaces  with  different  roughnesses  and 
ran  the  RBS  and  NR  spectra  through  the  standard  analysis.  As  we  reported 
in  reference  16,  the  RBS  distributions  were  consistent  with  assuming 
flat  iron  surfaces  with  small  amounts  of  carbon  in  the  outer  layers. 

Based  on  these  considerations,  we  neglect  surface  roughness  effects 
totally  in  our  calculations. 


15.  R.D.  Edge  and  U.  Bill,  MSurface  Topology  Using  Rutherford  Back- 
scattering,”  Nucl.  Inst.  $  Meth.,  168,  1980. 

16.  Andrus  Niiler  and  Robert  Birkmire,  "Measurement  of  Oxygen  and 
Nitrogen  Profiles  in  Steel,"  Nucl.  Inst.  $  Meth.,  149,  1978,  301. 
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There  is  a  maximum  depth  from  which  information  about  the  various 
elemental  concentrations  can  be  obtained.  This  depth,  referred  to  as 
analysis  depth  limit  in  subsequent  discussion,  depends  on  two  major 
factors;  the  deuteron  beam  energy  and  the  existence  of  interfering 
structures  in  the  RBS  or  NR  spectra.  The  beam  energy  dependence  arises 
because  the  particles  lose  energy  continuously  in  the  surface  layers,  and 
thus  a  depth  exists  which  is  large  enough  so  that  particles  with  a  given 
scattered  energy  cannot  emerge,  or  emerge  with  energies  too  low  to  be 
detected.  In  figure  3,  the  energy  cutoff  of  0.15  MeV  in  the  RBS  spec¬ 
trum  is  due  to  the  electronic  system  threshold.  Deuterons  detected  at 

2 

this  threshold  energy  were  scattered  from  a  depth  of  0.55  mg /cm  with 
scattering  from  deeper  layers  yielding  undetected,  lower  energy  deuter¬ 
ons.  The  fact  that  interfering  structures  can  affect  the  analysis 
depth  is  illustrated  by  the  way  the  nitrogen  peak  of  figure  3  runs  into 
the  oxygen  peak.  The  lowest  energy  at  which  the  nitrogen  does  not  in¬ 
terfere  with  the  oxygen  is  about  1.41  MeV.  The  protons  which  emerge 
from  the  sample  with  this  energy  were  produced  at  a  depth  of  about 
,  2 

0.53  mg/cm  . 


4.  THE  14N(d,p5)15N  CROSS  SECTION  MEASUREMENT 

The  most  commonly  used  reaction  for  nitrogen  analysis  is  the 
14N(d,a1)  12C. 17  However,  the  14^^^  reaction  cross  section  is 

about  four  times  higher  than  the  (d,a^)  and  since  the  reaction  gives 

the  highest  sensitivity  for  nitrogen  detection,  it  has  been  used  in 

this  experiment.  Although  the  (d,p  )  cross  section  has  been  deter- 

18**  19 

mined  at  energies  above  0.6  MeV  ,  the  data  is  sparse  below  that  energy 

and  not  adequate  for  ion  beam  analysis.  In  this  section,  the  measurement 
of  the  (d,Pj.)  excitation  function  between  0.32  and  1.45  MeV  deuteron  ener¬ 
gy  and  scattering  angle  of  160°  in  the  laboratory  will  be  described.  The 
work  in  this  section  has  been  published  in  reference  20. 

Two  accelerators  were  used  to  collect  the  data  in  this  experiment. 

The  energy  region  from  0.40  to  1.45  MeV  was  covered  in  0.05  MeV  steps 

TT.  J.W.  Mayer  and  E.  Rimini,  Ion  Beam  Handbook  for  Material  Analysis, 
Academic  Press,  New  York,  1977. 

18.  G.  Amsel,  and  D.  David,  "Microanalysis  of  Nitrogen  by  Means  of  the 
Direct  Observation  of  Nuclear  Reactions  Applications,"  Rev.  Phys. 

Appl .  4,  1969,  383. 

19.  A  Valek,  T.  Vertse,  B.  Schlenk  and  I.  Hunyadi,  "A  Study  of  the 

44N(d,p)4^N  Reaction  at  Low  Bombarding  Energies,"  Nucl.  Phys.,  A270, 
1976,  200.  5 

20.  A.  Niiler  and  R.  Birkmire,  "The  N(d,p5)  N  Cross  Section,  0.32-1.45 

MeV,"  Nucl.  Inst.  $  Meth.  168,  1980,  105. 
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at  the  University  of  Delaware  2  MeV  Van  de  Graaff  accelerator.  The  .872 

19 

MeV  resonance  in  the  F(p,ay)  reaction  was  used  to  calibrate  the  analyzer 
magnet.  Beam  slits  after  the  analyzing  magnet  were  set  so  as  to  give 
a  beam  with  energy  with  accuracy  better  than  ±  .002  MeV  throughout  the 
entire  operating  range.  The  BRL  Cockcroft-Walton  accelerator  was  used  to 
cover  the  energy  range  0.32  to  0.70  MeV  in  energy  increments  of  0.005  MeV. 

19 

The  .340  and  .669  MeV  resonances  in  the  F(p,cty)  reaction  were  used  to 
calibrate  the  analyzer  magnet  on  this  machine.  The  .006  MeV  width  of  the 
.669  MeV  resonance  and  long  term  drifts  in  the  system  limited  the  Cockcroft- 
Walton  energy  accuracy  to  no  better  than  ±  .005  MeV. 

Typically,  100  to  200  nA  of  deuterons  were  used  during  the  measure¬ 
ments.  Targets  were  self  supporting  foils  made  of  a  layered  structure 

2 

of  copper,  melamine  and  copper,  all  evaporated  onto  a  20  yg/cm  carbon 

foil.  The  total  target  thickness  was  estimated  to  be  less  than  150  yg/ 

2 

cm  .  Although  this  target  structure  gave  the  best  stability  under  beam 
heating  of  all  targets  tested,  the  targets  did  degrade  during  the  experi¬ 
ment.  Later  discussion  will  show  that  absolute  target  stability  was  not 
necessary.  Spectra  were  taken  only  at  160°  in  the  laboratory.  At  this 
angle,  it  was  possible  to  arrange  the  beam-target-detector  geometry  in 
such  a  way  that  the  incident  and  outgoing  particles  passed  through  only 
2 

one  30-40  yg/cm  copper  layer  on  their  way  to  and  from  the  melamine. 

The  incident  deuteron  energy  losses  in  the  copper  layer  were  at  most 
.01  MeV  and  all  reported  bombarding  energies  have  this  correction  factor 
applied.  In  addition,  the  melamine  layers  were  thin  enough  so  that 
neighboring  (d,p)  peaks  from  nitrogen,  oxygen  and  carbon  did  not  inter¬ 
fere  with  the  analysis. 


A  sample  of  the  spectra  taken  with  a  high  resolution  (12  keV  for 
5.5  MeV  alphas)  Silicon  Surface  Barrier  detector  is  shown  in  Figure  7. 
Fast  pulse  pile-up  rejection,  75  nsec  resolving  time,  was  used  in  order 
to  keep  pile-up  signals  from  copper  from  producing  additional  backgrounds 
under  the  (d,p  )  peak. 


Since  the  targets  did  not  remain  stable  under  beam  heating,  abso- 

14 

lute  cross  sections  for  the  N(d,ps)  reaction  could  be  obtained  only 
by  comparison  to  another  14N  cross  section.  The  unresolved  14(Nd,p1+2) 
peak  was  used.  The  (d,p^)  cross  section  is  given  by 


r(d,p5)  a(d,P1+2)  Y5/Yl+2 


(5) 


where  the  o' s  are  the  differential,  laboratory  cross  sections  and  the 
Y's  are  the  yields  of  the  p^  and  p^+2  groups.  In  addition,  neither 

the  solid  angle  or  beam  charge  measurements  are  needed  to  determine  the 
cross  sections  by  this  method. 
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Figure  7.  Spectrum  of  particles  produced  by  0.65  MeV  deuterons  at  160°  from  the  Cu-Melamine-Cu- 
Target . 


Values  for  the  reference  a,,  .  cross  section  were  obtained 

(d,Pl+2) 

18  21 

from  a  smooth  curve  drawn  through  the  available  data  of  Amsel  ,  Booth  , 
19 

and  Valek  .  Valek's  data  was  followed  below  .6  MeV  and  Booth's  and 

22 

Amsel 's  data  was  used  above  .6  MeV.  Although  Porto  has  performed 
measurements  at  energies  above  1  MeV,  it  could  not  be  read  from  the  graph 
with  enough  accuracy  for  present  purposes  and  was  thus  not  used.  Includ¬ 
ing  the  quoted  8%  uncertainty,  this  type  of  interpolation  procedure  in¬ 
creases  the  final  uncertainty  on  the  reference  cross-sections  to  about 
10%.  Over  most  of  the  energy  range,  this  10%  is  the  predominant  error  in 
the  measured  values.  However,  at  the  lower  energies,  the  counting  sta¬ 
tistics  become  more  important  so  that  the  worst  case  uncertainty  on  the 
cross  sections  is  about  35%. 

The  results  of  this  experiment  are  shown  in  Figure  8  along  with 
data  from  Amsel,  Valek,  and  Porto.  The  Van  de  Graaff  points  are,  in 
general,  from  a  single  data  run.  However,  since  the  Cockcroft -Walton 
energy  uncertainty  was  about  .005  MeV,  the  data  points  are  averaged 
over  each  .01  MeV  energy  interval.  It  can  be  seen  that  the  present 
data  generally  agree  within  the  uncertainty  limits  with  the  published 
data  of  Valek  and  Amsel.  The  fact  that  the  agreement  with  Porto's  data 
and  Valek's  point  at  .626  MeV  is  not  as  good  can  be  attributed  to  the 
(d,Pi+2)  normalization  data  that  we  used.  The  fact  that  all  the  data 

shown  are  not  at  160°  is  not  significant  since  the  angular  distributions 
from  Porto  at  1  MeV  and  up  seem  to  be  quite  flat  at  these  backward 
angles .  The  small  amount  of  structure  indicated  by  the  present  data 
near  0.6  MeV  has  not  been  observed  previously  although  Amsel 's  data 
does  flatten  out  just  before  reaching  .6  MeV.  These  new  data  improve 
the  definition  of  the  shape  of  the  (d,p5)  cross  section  below  0.7  MeV. 

5.  RESULTS  AND  DISCUSSION 


5.1  Concentrations 


Following  the  procedures  outlined  in  Section  3.1,  concentration 
profiles  were  obtained  for  twelve  nozzles.  The  concentrations  which 
gave  the  best  fits  to  the  experimental  spectra  are  shown  in  figures 
9(a)  through  20(a)  while  the  fits  are  shown  in  figures  9(b)  through 
20(b).  A  glance  at  figures  9(b)  -  20(b)  can  be  used  to  obtain  qualita¬ 
tive  information  about  the  amounts  of  C,N  and  0  on  the  samples.  The 
wider  peaks  indicate  greater  depths  of  penetration  as  can  be  seen  by 
comparing  figures  9(b)  (Noz  2)  and  17(b)  (Noz  10),  and  the  higher  yields 

21.  D.L.  Booth,  F.V.  Price,  D.  Roaf  and  G.L.  Salmon,  "The  Differential 

14  15  14  12 

Cross  Sections  for  the  Reactions  N(d,p)  N  and  N(d,a)  C  Between 
600  and  1000  keV, "  Proc.  Phys .  Soc.  71,  1959,  325. 

22.  V.  Gomes  Porto,  N.  Ueta,  R.A.  Douglas  and  0.  Sala,  "Deuteron  Induced 

Reactions  on  ^^N,"  Nucl .  Phys.  A136,  1969,  385. 


29 


=  0.65  MeV 


30 


the  best  estimate  of  the  correct  cross  section  and  shows  the  values  used  in  the  analysis. 
The  error  bars  indicate  both  counting  statistics  and  systematic  uncertainty  from  the 
reference  cross  sections.  Error  bars  were  omitted  from  some  points  for  clarity  but  their 
size  can  be  gauged  from  other  points  near-by. 
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Figure  9a.  Concentration  depth  profiles  for  nozzle  #2. 
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Figure  10a.  Concentration  depth  profiles  for  nozzle  #3. 
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Figure  11a.  Concentration  depth  profiles  for  nozzle  #4. 
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Lgure  12a. 


Concentration  depth  profiles  for  nozzle  #5. 
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Figure  12b.  Fit  to  the  RBS  and  NR  Spectra  from  nozzle  #5  with  concentration  profiles  shown  in 
figure  12(a). 
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Figure  13a.  Concentration  depth  profiles  for  nozzle  #6. 
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Figure  14a.  Concentration  depth  profiles  for  nozzle  7. 
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Figure  15a. 


Concentration  depth  profiles  for  nozzle  8. 
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Figure  15b.  Fit  to  the  RBS  and  NR  spectra  from  nozzle  #8  with  concentration  profiles  shown  in 
figure  15(a). 
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Figure  17a.  Concentration  depth  profiles  for  nozzle  #10. 
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Figure  17b.  Fit  to  the  RBS  and  NR  spectra  from  nozzle  #10  with  concentration  profiles  shown  in 
figure  17(a). 
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Figure  18a.  Concentration  depth  profiles  for  nozzle  #11. 
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Figure  18b.  Fit  to  the  RBS  and  NR  spectra  from  nozzle  #11  with  concentration  profiles  shown  in 
figure  18(a). 
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Figure  19a.  Concentration  depth  profiles  for  nozzle  #12. 


0.7 


51 


NOZZLE  12 


(pa)  ai3iA 


c 

•H 

G 

£ 

O 

■C 

m 

m 

<D 


4h 

O 

P4 

G 

O 


aJ 

+-> 

G 

<D 

U 

G 

O 

O 


CN 


N 

N 

O 

G 

e 

o 

tw 

aj 

-P 

O 

0 

P4 

to 

C* 

n3 

G 

aJ 

on 
eQ 
cc  , 


<D 

x: 


CTi 


O  <D 
P  P 
P 

+j  txo 
•H  *H 

tL,  4-1 


rQ 

CTi 


0 

P 

a 

•H 

PU 


52 


100 


80  - 


H 


r-1 


o 


60  - 


< 

Q£ 


NOZZLE  13 

-  CARBON 

-  OXYGEN 

- NITROGEN 

- IRON 

.  SULFUR 


LU 

U 


u 


20- 


ANALYSIS 

DEPTH 

LIMIT 


=fc=t 


n — ■ — i - 1 - r - 1 - r~ 

0.1  0.2  0.3  0.4  0.5  0.6 

DEPTH  (mg/cm2) 


Figure  20a.  Concentration  depth  profiles  for  nozzle  #13. 
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Figure  20b.  Fit  to  the  RBS  and  NR  spectra  from  nozzle  #13  with  concentration  profiles  shown  in 
figure  20(a). 


indicate  greater  concentrations.  However,  care  should  be  exercised  in 
the  latter  observation  because  the  yields  as  shown  are  relative  so  that 
only  the  NR  or  RBS  parts  for  a  given  nozzle  can  be  compared.  There  is 
no  common  scale  factor  between  spectra  from  one  nozzle  to  another.  The 
concentrations  shown  by  figures  9a  -  20a  can,  of  course,  be  compared 
with  one  another. 

There  are  some  common  features  which  are  observed  for  most  of 
the  nozzles.  Carbon,  nitrogen,  oxygen,  sulfur  and  barium  were  observed 
in  most  of  the  nozzles  surfaces.  The  carbon  and  oxygen  concentrations 
peak  at  the  surface  for  all  nozzles,  and  drop  off  to  relatively  lower 
values  by  the  time  the  analysis  depth  limit  is  reached.  The  nitrogen 
concentrations  generally  reach  their  maximum  at  a  small  distance  into 
the  surface  and  then  drop  off  to  negligible  values.  Sulfur  and  Barium 
are  found  only  on  the  outer  surface  layers. 

There  are  four  possible  sources  of  carbon  on  these  nozzles.  The 
first  is  from  the  carbon  used  to  mix  with  the  iron  to  form  steel  and 
should  amount  to  0.3  -  0.4%.  The  second  carbon  source  is  the  oil  which 
is  used  to  quench  the  steel  during  heat  treatment.  The  carbon  values 
from  the  surface  to  the  analysis  depth  limit  are  about  an  order  of 
magnitude  greater  than  in  the  bulk  and  are  due  to  the  quenching  oil.  The 
third  source  of  carbon  is  from  the  burning  propellant  gases  and  the  fourth 
possibility  is  from  vacuum  system  carbon  cracking  onto  the  outer  layers 
upon  the  surface  being  heated  by  the  deuteron  beam.  This  latter  source 
of  carbon  is  most  readily  observed  in  figure  9(a)  which  shows  a  100% 
carbon  peak  for  the  outermost  layer  on  nozzle  #2.  This  nozzle  was  used 
for  many  calibration  runs,  sometimes  in  relatively  poor  vacuums,  so  this 
surface  carbon  peak  is  not  surprising.  In  fact,  visual  observation  shows 
a  dark  spot  where  the  beam  hit  the  nozzle  surface.  Since  it  is  not  un¬ 
ambiguously  determinable  how  much  observed  carbon  is  due  to  heat  treat¬ 
ment  or  beam  heat  cracking,  the  contribution  from  the  propellant  gases 
cannot  be  determined  either.  A  more  careful  experiment  in  the  future 
using  better  vacuum  systems  and  defining  the  carbon  profile  of  a 
virgin,  unfired  surface  accurately  would  allow  a  better  determination 
of  carbon  content  changes  due  to  propellant  gases. 

Since  no  oxygen  is  allowed  in  the  base  steel,  there  are  only  two 
sources  for  it  in  the  nozzles.  The  first  is  simply  room  temperature 
oxidation  of  the  iron  at  the  surface  and  the  second  is  due  to  the  hot 
propellant  gases.  All  nozzles  show  a  maximum  value  for  oxygen  of  ~10% 
at  the  surface.  This  value  agrees  reasonably  well  with  oxygen  concen¬ 
tration  maxima  for  various  steel  samples  and  nozzles  which  have  not  been 
exposed  to  propellants.  Generally,  the  depth  of  the  oxygen  peak  is 

2 

less  than  .03  yg/cm  for  these  various  samples.  The  assumption  is  thus 

*  2 

made  that  all  oxygen  found  deeper  than  .03  yg/cm  is  due  to  the  burning 
propellants  while  it  would  not  be  possible  to  assign  a  precise  source 
to  the  oxygen  found  at  shallower  depths. 
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The  nitrogen  concentrations  generally  do  not  peak  at  the  surface, 

2 

but  at  a  distance  of  .03  to  .05  yg/cm  into  the  surface.  Also,  no 
nitrogen  is  found  on  unfired  nozzles  or  other  steel  samples.  This 
indicates  that  all  nitrogen  observed  is  due  to  the  hot  propellant  gases. 

The  reason  for  the  small  reduction  in  the  nitrogen  concentration  towards 
the  outer  surface  layers  cannot  be  explained  at  present. 

The  sulfur  and  barium  are  most  likely  due  to  the  BaSO^  and  KSO^ 

found  in  the  propellants  and  primers.  It  should  be  pointed  out  that 
part  of  what  is  reported  as  sulfur  could  in  fact  be  potassium  since  the 
RBS  with  deuterons  cannot  easily  distinguish  between  sulfur  and  potassium 
in  the  presence  of  the  much  larger  amounts  of  iron  and  no  convenient 
nuclear  reaction  peaks  exist  which  uniquely  identify  these  two  elements. 
Barium  concentrations  were  generally  in  the  few  tenths  percent  range  and 
are  not  shown  in  the  figures. 

The  concentrations  of  the  elements  C,N,0,S,  and  Ba,  as  determined 
by  applying  Eq.  4  to  the  analysis  depth  limit  are  given  in  Table  3. 

The  uncertainties-  on  these  table  values  are  approximately  15%  as  dis¬ 
cussed  in  Section  3.2.  From  figures  9(a)  through  20(a),  it  can  be  seen 
that  in  all  cases,  the  oxygen  and  carbon  extends  at  least  as  deeply  in¬ 
to  the  surface  as  the  analysis  depth  limit  but  in  many  cases,  the  nitro¬ 
gen  does  not  extend  so  far.  Since  all  oxygen  at  the  analysis  depth 
limit  and  beyond  would  be  expected  to  be  introduced  by  the  propellant 
exposure,  there  is  no  clear  way  to  determine  ambiguously  the  depth  of 
the  exposure  affected  layer.  In  order  to  do  so,  deuteron  beams  with 
higher  energies  should  be  used. 

5 . 2  Correlations . 

As  the  steel  surfaces  are  exposed  to  the  hot  propellant  gases,  many 
processes  can  contribute  to  the  production  of  the  oxide/nitride  layers. 

Some  of  these  are  absorption,  adsorption,  surface  reactions,  chemical 
reactions  and  diffusion.  Since  significant  amounts  of  oxides/nitrides 

are  found  at  least  as  deep  as  7500A  (several  hundred  atomic  layers)  a 
mechanism  such  as  preferential  diffusion  along  grain  boundaries  is  proba¬ 
bly  also  present.  The  effects  that  these  processes  have  are  complicated 
by  the  presence  of  the  shear  forces  of  the  gas  flow  which  act  simultaneous¬ 
ly.  Several  correlations  have  been  found  which  show  evidence  for  two  dis¬ 
tinctly  different  regimes  of  surface  erosion  caused  by  hot  gas  flows. 

Figure  21  shows  the  correlation  between  the  0+N  concentrations  and 
the  flame  temperatures  of  the  propellants.  The  four  points  at  each  of 
the  temperatures  represent  the  0+N  concentrations  on  the  nozzles  fired 
with  the  different  additive  conditions.  As  can  be  seen,  the  0+N  con¬ 
centrations  decrease  with  increasing  temperature  for  all  combinations. 

The  lines  are  drawn  to  indicate  the  strong  correlations  that  exist  for 
each  additive  type,  but  are  not  actual  fits.  Except  for  the  0+N  levels 


56 


TV  (°K) 


Figure  21.  Correlations  between  0+N  concentrations  and  the  propellant 
flame  temperatures . 
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Table  3.  Ion  Beam  Analysis  Results 
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are  driven  mainly  by  temperature  and  to  a  lesser  extent  by  the  additive 
type.  Sample  2  is  off  the  line  defined  by  10  and  6  because,  as  will  be 
discussed  later,  it  already  has  the  minimum  possible  0+N  level.  Figure 

22  shows  the  correlation  between  the  0+N  concentrations  and  the  measured 
mass  losses.  The  straight  line  curve  shown  fits  nine  of  the  data  points 
with  a  correlation  coefficient  of  0.90.  The  0+N  concentrations  decrease 
with  increasing  wear.  Figure  23  shows  the  correlation  between  the  pro¬ 
pellant  flame  temperatures  and  the  mass  losses.  The  straight  line  curve 
shown  fits  nine  of  the  data  points  with  a  correlation  coefficient  of 
0.91.  As  can  be  seen,  mass  loss  (or  wear)  increases  with  temperature. 
The  three  data  points  which  were  excluded  in  the  fits  in  figures  22  and 

23  were  #2,  6  and  8.  The  reason  for  this  exclusion  will  be  given  later. 
Figures  21-23  show  only  the  subsurface  0+N  concentrations. 

The  information  presented  in  figures  21-23  forms  the  basis  for  the 
following  model  for  erosion  from  this  set  of  nozzles.  All  exposures 
produce  oxygen  and  nitrogen  layers  on  and  in  the  surface  of  the  steel . 

23 

Recent  work  shows  tht  the  0+N  concentration  leve]  for  unfired  nozzles 
17  2 

is  about  1  X  10  at/cm  .  Three  nozzles  #2,  4  and  5  show  0+N  at  this 
"background"  level,  #6  is  barely  above  this  "background",  while  all  the 
others  show  significant  levels  above  "background".  Although  the  0+N 
layers  may  not  have  lower  melting  points,  they  certainly  have  less 
physical  strength  than  bare  steel  and  thus  are  more  easily  wiped  away 
by  the  hot,  flowing  gases.  In  fact,  as  can  be  seen,  the  hotter  expo¬ 
sures  wear  more  material  (see  fig.  23)  while  leaving  less  of  the  0+N 
layer  behind  (see  fig.  21).  Thus  it  appears  that  the  surface  is 
"softened"  and  simultaneously,  some  of  the  softened  layer  is  blown  away, 
the  amount  depending  on  the  exposure  temperature.  Nozzles  #2,  4  and  5 
appear  to  lose  all  of  the  0+N  layers  while  #6  is  on  the  borderline. 

The  fact  that  nozzle  #2  is  off  all  the  correlation  curves  (figs.  22 
and  23)  suggests  that  something  different  is  happening  during  exposure 
with  the  M2/no  additive  exposure.  Under  an  optical  microscope,  the 
nozzle  #2  surface  appears  smoother  and  more  polished  than  other  nozzles 
suggesting  a  melted  appearance.  As  mentioned,  this  surface  has  no  0+N 
levels  above  background  and  its  mass  loss  is  3-7  times  higher  than 
from  the  other  nozzles.  This  exposure  is  certainly  the  hottest  of  all 
the  twelve  in  this  experiment.  A  reasonable  conjecture  is  that  this 
exposure  is  hot  enough  so  that  the  melting  of  not  only  the  0+N  layers, 
but  also  of  the  base  steel  itself  can  occur.  Figures  22  and  23  show 
that  nozzle  #6  tends  toward  being  off  the  correlations  in  the  same  way 
as  nozzle  #2  but  to  a  lesser  extent.  It  is  possible  that  nozzle  #6  has 
the  second  hottest  exposure  of  all,  placing  it  in  a  region  where  steel 
melting  is  just  beginning  to  occur. 
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to  be  published. 
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Figure  22.  Correlation  between  the  0+N  concentrations  and  the  measured 
mass  losses. 
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Figure  23.  Correlation  between  the  propellant  flame  temperatures  and 
the  measured  mass  losses. 
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The  erosion  mechanism  indicated  by  the  lower  temperature  exposures 
is  the  formation  of  oxide,  nitride  and  perhaps  carbon  layers  by  physical, 
thermal  and  chemical  processes  followed  immediately  with  blow-off  by 
high  temperature  and  velocity  gas  flow.  We  label  this  mechanism  thermo¬ 
chemical  erosion.  The  higher  temperature  exposures  (nozzles  2  and  6) 
indicate  a  simpler  melt  and  wipe-off  mechanism. 

The  reason  that  the  three  data  points  from  nozzles  2,  6  and  8  were 
excluded  from  the  correlation  fits  is  that  they  all  deviated  from  the 
correlation  lines  much  more  than  the  average.  Since  the  correlation 
lines  fit  the  data  points  which  describe  the  thermo-chemical  erosion 
mechanism,  the  nozzle  2  and  6  data  would  not  be  fit  if  that  mechanism 
did  not  apply  to  them.  The  probable  reason  that  #8  deviates  so  much 
from  the  curves  is  that  it  was  fired  5  times  rather  than  4  but  the  first 
shot  was  made  with  no  blow-out  disc  in  place.  Such  a  significant  devia- 

23 

tion  from  normal  shots  has  been  found  in  a  recent  experiment  to  affect 
not  only  concentration  levels  but  subsequent  mass  losses  also. 

Figures  21  and  23  show  possible  evidence  that  the  use  of  wear  re¬ 
ducing  additives  do,  in  fact,  build  up  non-steel  surface  layers  which 
help  to  protect  the  samples  against  erosion.  Figure  23  shows  that, 
without  exception,  erosive  wear  is  reduced  whenever  additives  are  used. 
Figure  21  shows  that  with  but  two  exceptions,  the  N+0  concentrations  are 
greater  whenever  additives  are  used  but  those  two  exceptions  (10/12  and 
2/5)  are  well  within  experimental  uncertainties  of  each  other.  The  T^C^ 

wax  samples  show  the  greatest  amount  of  residue  on  their  surfaces  even 
though  the  erosion  from  those  samples  is  not  significantly  different 
from  the  Talc  Wax  or  PUF  samples.  This  fact  may  be  important  in  the 
choice  of  additive  under  specific  conditions  when  barrel  fouling  may 
prove  a  problem. 

5 . 3  Scanning  Electron  Microscope  Observations 

Figure  24  shows  the  SEM  pictures  for  the  four  nozzles  2,6,9  and  11. 
Figure  24a  is  for  nozzle  #2  and  shows  extensive  cracking  of  the  surface 
with  the  . 5y  cracks  running  in  jagged  patterns  suggesting  grain  bound¬ 
aries  in  austentite.  A  similar  severity  of  cracking  is  found  for  all 
the  nozzles  fired  with  M2  propellant.  Nozzle  shows  far  less  exten¬ 
sive  cracking  than  #2  but  about  the  same  as  #7,  8  and  9  which  were  also 
shot  with  M30  propellant.  Figure  24c  is  taken  from  nozzle  9.  Figure 
24d,  taken  from  nozzle  #11,  shows  no  cracks  whatsoever,  but  does  show 
the  machining  grooves  have  not  yet  been  wiped  away.  All  nozzles  fired 
with  Ml  propellant  show  features  very  similar  to  those  seen  in  this 
picture  of  #11.  The  extent  of  surface  cracks  does  not  correlate  with 
the  amount  of  wear  but  does  correlate  to  the  exposure  temperature.  The 
cracking  is  thus  most  likely  due  to  the  rapid  thermal  expansion  and  con¬ 
traction  cycling  experienced  by  the  surface. 
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Figure  24. 


Scanning  electron  microscope  photographs  of  the  surfaces  of 
nozzles  a)  #2,  b)  #6,  c)  #9,  and  d)  #11. 
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6.  CONCLUSIONS 


It  has  been  shown  that  a  method  now  exists  which  is  capable  of 
detailed  characterization  of  atomic  scale  surface  composition  changes 
produced  by  the  exposure  of  steel  samples  to  propellant  combustion 
products.  This  ion  beam  analysis  method  has  depth  resolution  capability 
of  about  0.05  pm,  detectability  limits  in  the  order  of  0.5  atomic  percent 
and  an  absolute  accuracy  of  about  1  atomic  percent  for  concentrations  of 
oxygen,  nitrogen  and  carbon.  In  addition,  this  method  is  non-destructive, 
allowing  a  series  of  experiments  to  be  conducted  on  a  given  sample. 

Most  importantly,  it  has  been  shown  that  erosion  mechanisms  may  now  be 
identified  from  the  results  of  such  surface  analyses. 

The  present  experiment  has  shown  that  the  effects  of  propellant 
burns  on  the  surface  compositions  of  gun  steels  depend  strongly  on 
propellant  type  (i.e.  propellant  flame  temperature)  as  well  as  the  type 
of  wear  reducing  additive  used.  Strong  correlations  have  been  found  for 
the  total  oxygen  plus  nitrogen  concentration  levels  and  the  propellant 
flame  temperatures  as  well  as  the  amount  of  erosive  wear  observed.  The 
flame  temperatures  and  measured  wear  also  correlate  strongly.  The 
analysis  of  the  results  indicate  that  there  are  at  least  two  regimes 
of  propellant  gas  erosion.  The  first  regime,  at  lower  temperatures, 
is  characterized  by  the  production  of  oxide/nitride  surface  layers  and 
the  immediate  removal  of  part  of  the  produced  layer.  The  extent  of 
the  production  and/or  removal  depends  on  the  propellant  flame  tempera¬ 
ture  such  that  cooler  propellants  leave  higher  oxygen/nitrogen  levels 
than  the  hotter  ones.  This  mechanism  is  most  likely  thermo- chemical 
in  nature.  The  second  regime,  at  higher  temperatures,  is  characterized 
by  the  complete  absence  of  oxides  or  nitrides  above  background  levels 
as  well  as  a  surface  that  is  smooth  and  well  polished  in  appearance. 

This  mechanism  is  most  likely  surface  melting  and  immediate  wiping 
off  of  the  melted  material  by  the  gas  flow  forces. 

There  are  several  open  questions  left  by  this  work.  Since  the  ion 
beam  analysis  is  sensitive  only  to  the  atomic  concentrations,  no  detailed 
chemical  information  can  be  obtained  from  the  data.  Other  methods  such 
as  x-ray  photoelectron  spectroscopy  (XPS)  or  secondary  ion  mass  spectro¬ 
scopy  (SIMS)  can  be  used  to  definitively  determine  the  chemical  forms 
that  the  observed  oxygen  and  nitrogen  are  in.  Such  information,  coupled 
with  the  present  results,  could  then  be  used  to  better  define  the  thermal 
properties  of  the  produced  layers .  Any  observed  differences  in  the 
chemical  compositions  of  the  surface  layers  would  lead  to  improved 
understanding  of  the  erosive  effects  of  individual  propellants.  A 
more  precise  measurement  of  changes  in  the  carbon  concentrations  would 
also  be  important.  The  extent  of  initial  and  exposure  related  carburi¬ 
zation  is  highly  relevant  to  the  strength  of  the  steel  surface  and  its 
ability  to  resist  erosion. 
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